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Abstract 

In this paper we study a new class of shell models, defined in terms 
of two complex dynamical variables per shell, transporting positive and 
negative helicity respectively. 

The dynamical equations are derived from a decomposition into helical 
modes of the velocity Fourier components of Navier-Stokes equations (F. 
Waleffe, Phys. Fluids A 4, 350 (1992)). 

This decomposition leads to four different types of shell models, ac- 
cording to the possible non-equivalent combinations of helicities of the 
three interacting modes in each triad. 

Free parameters are fixed by imposing the conservation of energy and 
of a "generalized helicity" in the inviscid and unforced limit. 

For a = 1 this non-positive invariant looks exactly like helicity in 
the Fourier-helical decomposition of the Navier-Stokes equations. 
Long numerical integrations are performed, allowing the computation 
of the scaling exponents of the velocity increments and energy fiux mo- 
ments. 

The dependence of the models on the generalized helicity parameter a 
and on the scale parameter A is also studied. 

PDEs are finally derived in the limit when the ratio between shells goes 
to one. 
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1 Introduction. 



The most intriguing problem in 3 dimensional turbulence is related to the un- 
derstanding of the dynamical mechanism triggering and supporting the energy 
cascade from large to small scales. According to the celebrated Kolmogorov 
theory (K41), energy should be transferred downwards in scales following a 
self-similar and homogeneous process. This idea, plus the assumption of local 
isotropy and small scales universality, led Kolmogorov to formulate a precise 
prediction on the scaling properties of the increments of turbulent velocity 
fields, 6v{l) ~ \v{x + I) — v{x)\, at distances /. Namely, following Kolmogorov 
and denoting by e(/) the energy transfer through scale /, we have: 

S,{1) ^ {{6v{l)y) = C,{{e{l)y/')F/' ~ /«^) (1) 

where Cp are constants and the scale / is supposed to be in the inertial range, 
i.e. much smaller than the integral scale and much larger than the viscous 
dissipation cutoff. From (0) we deduce that depending on the scaling properties 
of the energy transfer, we can observe different power law behaviors for the 
structure functions Sp{l). Kolmogorov made the simplest assumption by taking 
e(/) to be constant and deducing that: 

C(p)=p/3,Vp (2) 

While from a qualitative point of view the Kolmogorov intuition was a true 
breakthrough in the understanding of turbulence, his theory lacks quantita- 
tive agreement with experiments. In particular, there are many experimental 
and numerical 0, |^, |^ results telling us that energy is transferred intermit- 
tently, leading to non-trivial scaling corrections to the "p over 3" Kolmogorov 
prediction for the ({p) exponents. 

More and more efforts have been devoted in trying to understand and to model 
with good accuracy the energy transfer dynamics. 

Beside analytical and direct numerical approach, there are two other possible 
choices: either building simple random processes for the chaotic energy transfer 
among different scales |^, |^, ^ or studying a dynamical deterministic model. 
In the last 20 years, among the possible dynamical modelizations of fully devel- 
oped energy transfer, those obtained by using shell models have been particu- 
larly successfully. These models concentrate only on the dynamical interactions 
among degrees of freedom at different scales, neglecting completely their spatial 
location. In this way a one-dimensional chain of interacting Fourier modes is 
constructed. The same non-linear structure of Navier-Stokes eqs. is retained, 
but their three-dimensional vectorial nature is completely lost. The simplifi- 
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cations are so strong that only a posteriori one can say whether the model is 
interesting and reliable or not. 

For the most famous among all shell models, the Gledzer-Ohkitani-Yamada 
(GOY) model, are now available a wide number of works and contributes (0- 
||15|| ). The intermittent and chaotic dynamics of the GOY model is remarkably 
close (for a suitable choice of the free parameters) to what is found for the 
Navier-Stokes eqs. 

The most striking property of the GOY model is that, for some particular 
choices of the free parameters, it has an intermittent and chaotic dynamics 
quantitatively close to what is found for the Navier-Stokes equations. In par- 
ticular, the set of scaling exponents ({p) coincides (for a suitable choice of 
parameters) with that measured in true turbulent flows. 

Recently, it has been pointed out that the GOY capability of having a chaotic 
and intermittent behavior with scaling exponents very close to the real ones, 
is strictly related to the conservation in the inviscid, unforced limit of two 
quadratic quantities. The first quantity is the energy, while the second is, 
roughly speaking, the equivalent of helicity in 3D turbulence [0]. 
Both the GOY-helicity ||l5l and the helicity in the N-S equations [T^, [1^ have 
been suggested to play an important role in triggering the intermittent nature 
of the energy cascade. 

This new emphasis on the role played by the helicity motivated us to study a 
new class of shell models. These models can be seen as a generalization of the 
GOY model such that their helical structure is closer to the corresponding N-S 
case (indeed, one of them is nothing but the old GOY case). 
We study shell models with two complex variables per shell, derived on the 
basis of an helical formulation of the N-S equations in the Fourier space. 
In this way, it is possible to obtain a second non-positive defined invariant 
closer to the definition of helicity in the N-S equations. Our aim consists in 
trying to understand the importance of both (competing?) transfer of energy 
and helicity in N-S eqs. by examining the non trivial dynamics shown by this 
new class of helical-shell models. 

The outline of the paper is the following: in section 2 we review the GOY model; 
in section 3 we introduce the new class of helical-shell models; in section 4 the 
basic triad interactions within three contiguous shells are studied. Section 5 
contains the results of our numerical simulations; in section 6 PDEs for the 
continuum limit (ratio between shells that goes to one) is derived; conclusions 
follow in section 7. 
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2 The GOY model. 



The GOY model can be seen as a severe truncation of the Navier-Stokes equa- 
tions: it retains only one complex mode representative of all Fourier 
modes in the shell of wave numbers k between kn = ^qA" and fcn+i, A being an 
arbitrary scale parameter (A > 1), usually taken equal to 2. 
The dynamics is governed by the following set of complex coupled ODEs where 
only couplings with the nearest and next nearest shells are kept: 



— ikn (^aU^^iUj^_^_2 + hu^j^iU^_i + Ctt„_^n„„2) ^k^Jln + Sn,nof^ (3) 

where u is the viscosity, / the external forcing acting on a large scale uq and 
a,b,c are three free parameters. By adjusting the time scale we can always fix 
a = 1; the possible choices for b,c are restricted by imposing the conservation 
of two quadratic quantities in the inviscid and unforced limit: 

W^l,2 = E4,2l«nr, (4) 
n 

where Zi^2 are the solutions of the quadratic equation: 

cX'^z^ + b\z + a = 0. (5) 

In order to stay as close as possible to the N-S equations, we require that one 
of the two conserved quantities is the energy, i.e. zi = 1: 

Wi = E = J2K\^- (6) 

n 

If we rewrite: 



Xb = -e 



2._. 1 (7) 



we are left with only two (A and e) of the original four free parameters. The 
second quadratic invariant is: 



W2 = H = Y.{e- 1) 



The characteristics of this second invariant change by changing e: when e < 1 
it is not positive-defined (as helicity in 3D), while if e > 1 it is positive defined 
(as enstrophy in 2D). Expression (§) can be rewritten as: 
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^« = Ex(e)X^^''VnP, (9) 

n 

where x{^) = sign{e — 1) and the a parameter is related to e and A by: 

|e-l| = A-". (10) 

Our interest here is to consider how well the dynamics of a 3D turbulent flow is 
reproduced by the model: only e in the range < e < 1 will be taken, in order to 
have a non-positive defined second invariant. Indeed, for a = 1 our "generalized 
helicity" Ha = Z]n(~l)"^n I'^^nP physical dimensions coinciding with the 3D 
Navier-Stokes helicity. The two free parameters of the model can be taken to 
be A (the ratio between adjacent shells) and a. The two coefficients b and c 
can be rewritten as: 

b = x-"^-^ - \-\ , 

Such a class of models has a highly non trivial dynamical behavior. Intermit- 
tency of the energy transfer and multifractal nature of energy dissipation have 



been studied in R 10, 11, 12, 13 



It turns out that the values of the C(p)s are not universal, depending on the 



choice of e and A 0, |T4|. Nevertheless, Kadanoff et al. [|TJ] verified that the 

scaling exponents are invariant along the curve in the (e. A) plane where both 

energy and helicity are conserved, i.e. the curve at a(e. A) = 1. 

This suggested that the second invariant plays a crucial role in the model 

dynamics. 



More recently Biferale and Kerr [|T5| have attributed to the helicity the role of 
triggering the intermittent cascade of the energy from large to small scales. 
These considerations together with the observation that this "GOY-helicity" 
is only partially consistent with the N-S helicity (i.e. it presents an asymmetry 
between odd and even shells that does not have any counterpart in physical 
flows) persuaded us to study a modified shell model [jl5l, with two complex 



variables in each shell, carrying helicity of opposite sign, in order to obtain a 
second invariant closer to the N-S helicity. In the following section we introduce 
this new class of shell models, whose nonlinear interactions are constructed on 
the basis of a helical decomposition of the N-S equations in the Fourier space. 

3 The helical-shell models. 

In order to introduce two helical variables per shell we refer to the velocity field 
in N-S equations, expanded in a Fourier series |TH]. The velocity vector can be 
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represented in terms of its projection on an orthogonal basis formed by k, h+ 
and h_. The two basis vectors h_|_ and h_ can be chosen to be the eigenmodes 
of the curl operator: 

ik X hg = skhs, (12) 

where s = ±1. 

This corresponds to an expansion of the velocity vector into helical modes: 



u(x) = ^u(k)exp(ik-x) = ^[M+(k)h+ + M"(k)h_] exp(ik ■ x). (13) 

k k 

The real flow velocity corresponding to the plus (minus) mode rotates clockwise 
(counterclockwise) as one moves in the direction of k, thereby forming a left- 
handed (right-handed) helix; the vorticity vector of such a flow is parallel 
(antiparallel) to the velocity and the helicity is maximum (minimum). 
The kinetic energy and helicity are given by: 



E = Eki5;(k) = Ek |u(k) ■ u*(k) = Ek(|u+(k)P + |u-(k)n, 

H = Ek^(k) = Ek|u(k) ■ ^*(k) = EkM|u+(k)p - |u-(k)| 



(14) 



Plugging eq.(|T3D into the N-S equations yields to the dynamical evolution for 
the complex amplitudes u'^''(k,t) {sk = ±1) 



|«^Hk) + z/fc V^k) = E E 9KpA^pP-s,q){u^^{py^{q)r. (15) 

k+p+q=0Sp,S5 

The geometric factor 5'k.p,q = — |(hsi, A hgp ■ hg^)* can be developed and fac- 
torized: 

Skk + Spp + Sgq 
g = r , (16) 

P 



where r is a function of the triad shape only ||18[| . 

Eight different types of interaction between three modes M'^'=(k), u^p{p), u'^«(q) 
with |k| < IpI < |q| are allowed according to the value of the triplet (s^, Sp, Sg) 
= (±1,±1,±1): among them, only four are independent, the coefficients of 
the interaction with reversed helicities (— s^, —Sp, —Sq) being identical to those 

with {Sk,Sp,Sq) 



1. (sfc, Sp, Sg) = (+, -, +) or (-, +, -) 

2. (sfc, Sp, Sq) = (+, -, -) or (-, +, +) 
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3. (sfc, Sp, Sq) = (+, +, -) or (-, -, +), 



4. (sfc, Sp, Sq) = (+, +, +) or (-, -). 

The equations corresponding to the single interaction (s^, Sp, Sg) have the form 
(omitting viscosity and forcing): 

ii"' = r{spp - Sqg)^^^±^^^±^(u"PM"«)*, 
u^p = r{sqq - Skk)^i^^^^^f±^{u'm'^y, (17) 
= r{skk - Spp) "'^+'pP+'^i (u'^u'''y . 

Each interaction independently conserves both energy and helicity on a single 
triad. 

The dynamical system (0) composed by a single triad can be considered as 
the basic brick of the semi-infinite chain leading to the transfer of energy in 
turbulent flow. 

By studying its stability properties it is possible to understand how energy and 
helicity are transferred among different wave-vectors belonging to the same 
triad. 

Following |[T^, we distinguish two different kinds of dynamics: for the cases 
corresponding to the choices 1 and 3 of the triad helicities, the unstable wave- 
vector is the smallest-one, while for the cases 2 and 4 the unstable wave- vector 
is the medium-one. This very simple analysis suggests that by linking together 
a series of triads we should have a forward energy transfer for cases 1 and 3 
and both forward and backward (competing) energy transfers for cases 2 and 
4. 

In a turbulent flow the direction of energy transfer is dynamically controlled 
by the triple correlation < 'u'^'=(k)M*p(p)'u*«(q) >. It is reasonable to argue that 
the statistical properties of < u''^''{k.)u^p{p)u''^''{q) > are such that the overall 
direction in energy transfer coincides with the simplified behavior inferred from 
the stability study of the single triad {instability assumption in [|^). For 
instance, it is easy to estimate, by using the instability assumption, what would 
be the net energy transfer in the above four cases if the energy spectrum had 
the Kolmogorov scaling, E{k) = k~^^^ p8[]: 



1 and 3: direct energy cascade from large to small scales, 

4: reverse energy cascade from small to large scales, 

2: direct (reverse) energy cascade for local (nonlocal) triads. 
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The helical decomposition of the N-S eqs. suggested us the opportunity of 
defining a different GOY-like shell model for each one of the above four classes. 
In each shell we will have two complex dynamical variables and u~, trans- 
porting positive and negative helicity respectively: 

(18) 

where j = 1, 4 labels the four different models and the helicity indeces in the 
non linear interactions arc easily found for each of the four cases (see table(l)). 
The coefficients aj,bj,Cj are determined imposing, as usually, the energy con- 
servation: 

^E^^(E(wr+K-r))=o. (19) 

that leads to the same relation for all models: 

aj + bjX + CjX^ = (20) 
By imposing also the conservation of the generalized helicity: 

|^^«-|i:oi<p-Ki^)=o (21) 

we obtain different relations for the four models: 



1. ai - A°+i6i + A2("+i)ci = 0, 

2. 02 - A°+i62 - A2(°+i)c2 = 0, 

3. 03 + - A2("+i)c3 = 0, 

4. a4 + A°+^64 + A2("+^)c4 = 0. 



Fixing ttj = 1 one then finds the expressions for the coefficients bj and Cj in 
terms of the parameters A and a (see table(2)). 

Let us remark two important facts. First, model 1 is nothing but two masked 
and uncorrelated versions of the original GOY model, with dynamical variables 
{uf , ■■■) and (u^ , ■■■) respectively; rewriting the coefficients bi 

and Cl in terms of the usual parameters A and e one can easily recover the 
standard GOY model coefficients. Second, also model 4 is formed by two 
independent sets of variables {uf, u^, u^, ■■■) and {ui, U2, u'^, ■■■), each of them 
conserving separately a positive-definite quantity similar to enstrophy in 2D. 
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Thus, model 4 is equivalent to two uncorrelated GOY models for 2D turbulence 

EH- 

The fact that the model 1 is formed by two uncorrelated GOY models is clearly 
due to our choice of taking only first and second-neighbor interactions. Model 
4, on the other hand, will always be the sum of two separated models for any 
choice of the interacting modes composing the triads. 

In the following, we will refer to the properties of model 1 intending corre- 
sponding properties of the GOY model. Model 4 will be studied only for 
completeness. 



4 One-triad systems. 

Following the instability assumption |TB| that connects the single triad dynam- 
ics with the global statistical behavior of a multi-triads flow, we repeat the 
analog stability study for the three shells, single triad system. 
By isolating three shells of waves numbers ki,k2,ks, we can inspect their dy- 
namical properties as determined by their mutual interactions. 
For the positive-helicity modes we have: 

u+ = ik^iul^u^"')*, 

ut = tk2b,{u'^'^ul'y, (22) 
iit = iksCjiuful')*, 

where j = 1,...,4 for the four models. An analogous set of equations holds 
for the negative-helicity modes, changing the sign of the helicity index of all 
variables. 

This system conserves both energy and helicity. 

The corresponding equations for the energies are the following: 

Ei=A, 

E2 = b,XA, (23) 
E3 = c,\^A, 

where A = 2A;i5>[(m^^m^^m|) + (% ^^m^'^mD]. 



As found in [|18|, we know that the unstable mode is: 

• the smallest mode for interactions 1 and 3, 

• the medium mode for interactions 2 and 4. 
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In order to have a deeper understanding of the energy transfer dynamics, we 
have performed several integrations of eqs.(|^), using the parameters values 
X = 2, a = 1, ki = and different initial conditions. 
This analysis, performed on all four models, gives the following results: 

• model 1: mode 1 gives energy equally to mode 2 and mode 3, 

• model 2: mode 2 gives more energy to mode 3 and less to mode 1, 

• model 3: mode 1 gives more energy to mode 2 and less to mode 3, 

• model 4: mode 2 gives more energy to mode 1 and less to mode 3. 

These energy exchanges are summarized in fig.(l). 



Behaviors 1 and 4 have already been noticed by Ditlevsen and Mogensen ^1 
for the 3D and 2D GOY model respectively. 

It is also interesting to investigate how these properties are modified when 
varying the a parameter in the models. Considering that in eq. (^3]) the sum 
of the three right-hand sides must be zero and normalizing to one the energy 
rate on the unstable shell, one can evaluate how the energy sharing between 
the other shells is affected by changing a (see fig. (2)): 

• In model 1 there is a clear dependence on a. As this parameter grows 
up, more and more energy is captured by mode 2. For a > 1, the energy 
gained by mode 2 become greater than the energy gained by mode 3, 
leading to a more local energy transfer; 

• model 3 is remarkably independent of a, as we shall see in the following. 
This fact has very important consequences for the intermittent dynamics 
of the complete shell model; 

• model 2 has a trend analogous to that of model 1, but with more drastic 
consequences: at a ~ 1.27 the mode that receives most of the energy 
from the unstable mode becomes the first instead of the third-one. This 
would suggest a change in the direction of the flux: from downward to 
upward; 

• in model 4 the mode that receives most of the energy remains mode 1, 
for all values of a, with a consequent reverse energy flux in all cases (as 
it must be, being model 4 a couple of 2D GOY models). 
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What emerges from this analysis is that the behavior of models seems to depend 
on the choice of the free parameter a, sometimes with strong consequences (as 
the direction of the flux in model 2). The only remarkable exception is the 
very low dependence of model 3. 

Concerning helicity, we can consider the following equations: 

H,{a) = B, 

H2{a) = r/2&iA"+i5, (24) 

where B = 2ki'^^'^[{ul^U2^Ui ) — {u^^^U2^^u{)] and 772 and r]s depend on the 
particular model considered (see table(3)). 

By performing the same analysis done for the energy evolution, one can con- 
clude that helicity is transferred in different ways, as depicted in fig. (3). 
Being helicity a non-positive defined quantity, forward (backward) transfer 
of positive (negative) helicity is equivalent to backward (forward) transfer of 
negative (positive) helicity. In view of this trivial remark, arrows in fig. (3) have 
only a visual value: indicating how helicity (with its own sign) is redistributed 
among shells. 

Let us notice that models 1 and 3 show a very different pattern in the helicity 
exchange among shells. This can be the explanation of the very different scaling 
properties shown by the two models when varying a (see next section). 
For example, we could argue that the dramatic dependence of the energy ex- 
change on the a parameter in models 1 and 2, together with the well defined 
direction of the helicity transfer, can somehow enhance the role played by the 
second invariant with respect to the other two models. 

In the next section we will check these arguments by performing a numerical 
study of eq.(|T|). 

5 Numerical analysis (model 3). 

In this section we will concentrate on the study of the statistical properties of 
model 3 compared with the already known results for the GOY model (model 

1)- 

Model 3, at difference from model 2, shows for any value of a a forward energy 
transfer. 

We integrated eqs.(|l8|) for model 3 using the standard parameters a = 1, A = 2, 
ko = 2~^, = 5(1+2)10"^, no = 1, u = 10"'' and a total number of shells equal 
to = 22 and = 26. In the numerical integration we used a fourth-order 
Runge-Kutta method, with a time step varying between dt = 10~^ (for the 
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simulations with 22 shells) and 10~^ (for the cases with 26 shells). Most of the 
results presented here are for the case N = 22, with a number of iterations of 
the order of hundred millions, which correspond roughly to several thousands of 
eddy turnover times at the integral scale. Stationarity is checked by monitoring 
the total energy evolution. 

The quantities we have looked at are the structure functions: 

Spin) = {\unn, (25) 

where Un = y^jn+p + 

In fig. (4) we show \og2[Sp{n)] as functions of log2(A;n) for p = I, 8. There is 
a well defined inertial range where the structure functions follow a power law: 

Spin) ~ (26) 

Let us notice that in this model, at variance from the GOY model, there are 
not the period three oscillations superposed on the power-law scaling. In this 
case a linear least-square fit allows to compute the scaling exponents, C{p), 
with an uncertainty smaller than in the GOY model. 

Nevertheless, in order to have a better estimate of the C{p)s one can study the 
moments of a particular third-order quantity, the mean energy flux through 
the ?7,th shell: 



n(n) = 2fc„S5[(<V2<+iWn) + &j-(<+i«X'-i) + ^(<Vi<«^ 



As pointed out by Pisarenko et al. [13], one can write for this quantity the 
equivalent of the Kolmogorov's four-fifth law, expressing the balance between 
energy input and energy dissipation in the system. 
Considering the energy variation over the first n shells: 

J n n 

^[E ((!<!') + = -2^ E ^^((l«;t.P) + (fcP)) + 

m=l m=l 

+n(n) + 2g?[(/x;) + (/"<*)]> (28) 

and assuming a statistical steady state, in the limit of vanishing viscosity we 
are left with an inertial range in which n(n) is constant: 

n(n) ~ const. (29) 
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In analyzing the scaling properties of all our results we have always used Ex- 
tended Self Similarity (ESS). ESS consists in plotting one structure function 
versus another (5*3 (n), for example). ESS turned out to improve the precision 
with which scaling exponents can be measured in true turbulent flows p, Q 
and in shell models . 

We have applied ESS analysis to two kinds of fit: 

1. log^lSpin)] vs log2[S3(n)], 

2. log2 [Sp(n)] vs log2[S3(ra)], 
where 

Ep(n) = (|n(n)/A;„,r/3). (30) 

In all our simulations we have found the two sets of exponents coinciding within 
the numerical and statistical errors. 

Fig. (5) shows the C(p)s of model 3, compared with those of the GOY model 
with the same parameters a = 1 and A = 2 (corresponding to the classical 
choice which conserves the analog of the 3D helicity). The scaling is nearly 
the same. Indeed, we argued in the previous section that both models have 
a forward energy flux; what turned out to be different was the exchange of 
helicity among shells, together with a different sensitivity on the parameter a 
connected to this second invariant. 

What occurs is a strong similarity for a = 1; nevertheless we expect a different 
behavior when this parameter is allowed to vary. 

5.1 The a dependence 

For the a dependence of the models we have explored two other different values: 
a = 0.5 and a = 2 (keeping fix A = 2). 

It is well known |]12|, that the GOY model shows a strong dependence of 

its statistical properties on the a value. For example, if a < 2/3 the dynamics 
is attracted toward a fixed point with Kolmogorov scaling, ({p) = p/3. For 
a > 1 intermittency become more important than what is usually measured in 
turbulent flows [Q. 



On the other hand, the statistical properties of model 3 turn out to be robust 
under changes of the a parameter. 

In fig. (6) we show the (a{p) exponents for the GOY model and model 3 at 
a = 0.5,1,2. Clearly, there is an evident dependence of the C(p)s on a for 
the GOY case while for the model 3 the different exponents coincide within 
numerical errors. 
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This behavior is in perfect agreement with the phenomenological speculations 
argued in the previous section from the study of the single-triad system. 
The robustness of model 3, with respect to variation in a, gives to this model 
an important role among the possible shell models of turbulence. 



5.2 The A dependence 

Concerning the dependence on the scale parameter A, we have performed an 
exploratory study by fixing a = 1 and taking A = 1.5 and A = 2.5 in model 3. 
For the case with A = 2.5 the ({p) exponents are still stuck to the previous 
values at A = 2. On the other hand, for the case A = 1.5 we found a week 
discrepancy comparable with the one found in for the GOY model. The 



issue of what happens in the limit A ^ 1 (the so-called continuum limit) is one 
of the most intriguing problems that must be analyzed in both cases, GOY 
model and model 3 (see next section). 

In fig. (7) we show three sets of exponents obtained for different choices of A for 
both model 3 and the GOY model. 

The origin and significancy of the weak spreading in the values of C{p)s is far 
from being understood. By changing A, one changes the ratio between adjacent 
shells and therefore how viscous and inertial ranges match together. This non 
trivial matching maybe interferes also with the determination of the scaling 
exponents . 



The A dependence in all this shell models is however a very important open 
question due to the obvious interest in having a PDE describing the contin- 
uum limit (A 1) of the energy transfer. In the following section we derive 
the equation for the continuum limit of all four models and we present some 
proposal for further investigations. 



6 The continuum limit 

As anticipated in the previous section one of the most interesting and still 
unexplored aspect of GOY-like shell models is their dynamics in the continuum 
limit 

For continuum limit we intend the limit when the separation between shells 
goes to one, i.e. A ^ 1: 

kn+i = Xkn ~ (1 + 5)kn, (31) 
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where we have defined A = exp{6) ~ 1 + 5 + 0(5^). In the limit (|3T| ) we can 
expand the Un set as follows: 

u{K+m) = u{kn) + 6mkndku{kn) + 0(5^). (32) 

Taking into account the equivalent expansions for the a, 6, c coefficients ( see 
table 2) and after some simple algebra one realizes that all the three models 
(model 1,2 and 3) lead to the same expression in the continuum, namely: 

dtu^{k) = ik (4:ku~dkU'^ + 2ku'^dkU~ + {2 + a)u'^u~ — au~u~^ + 

-ueu+ + f{k), (33) 

the corresponding equation for the u~s is obtainable from (|33|) by changing all 
helicity indeces. 

This PDE worths a deeper study for many reasons. 

First, let us notice that the continuum limit is highly non reversible, i. e. 
trying to come back to a logarithmically-equispaced shell structure one does 
not recover the original equations (|TB|). Second, the continuum model shows 
an unexpected universality: it is the limit of three models which have very 
different behaviors at A > 1. Third, even in the continuum there are two 
conserved quantities (in the unforced and inviscid limit) corresponding to the 
continuum analogous of energy and generalized helicity: 

E = f ^(|n+p + \u-\'), H^ = I ^k"{\u+\' - \u~\'), (34) 



k'' ' ' ' ^' " J k 

where the, apparently unusual, dk/k integration step comes from the original 
logarithmically- equispaced shell structure. Let us remark that the most inter- 
esting difference between the continuum expression (^) and the analog for the 
old GOY model is that now in ( P^D also helicity conservation is well defined. 
This was not the case for the continuum GOY model. The apparent paradox 
(model 1 is formed by two uncorrelated GOY models when A > 1) is easily 
solved by noticing that in the continuum limit shells collapse in such a way 
that the original ordering is destructed. This limiting procedure introduces a 
coupling between the two sub-models. 

This drastic difference with the GOY continuum case suggests the possibility 
that this new set of PDEs has a much richer dynamics than the corresponding 
GOY PDEs. In that case, indeed, is quite easy to realize that PDEs are 
integrable along the characteristics |Q (at least for the case of real variables). 



The solutions have a burst-like shape with a Kolmogorov scaling, reaching 
infinite k at finite time (for zero viscosity). 
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Whether eqs. ( |33D are more interesting or not is still an open question. 

As for the continuum limit of model 4 we need to go the the second order in 

the 6 expansion and we obtain: 

dtu+{k) = tk{16k\dkU+f + 2k\+dlu+ + 12{a + 2)ku+dkU+ + 

+ {a + 2)\u+yy -iyk\+ + f{k). (35) 

In this case the continuum limit is much more similar to what should be the 
continuum limit of a shell model describing 2D turbulence. The only two 
conserved quantities are both positive-definite and coincide with energy and 
with a generalized enstrophy Qa- A much more detailed investigation of both 
models is postponed to a forthcoming study. 

7 Conclusions 

In this paper we have performed a detailed investigation of a new class of 
helical-shell models. From the helical-Fourier decomposition of Navier-Stokes 
eqs. we have extracted four non-equivalent types of shell models having two 
inviscid quadratic invariants similar to the conserved Navier-Stokes quantities. 
Two of these four models coincide with the 3D and 2D versions of one of the 
most interesting historical shell models: the GOY model. On the other hand, 
the other two (model 2 and model 3 in the text) show different and peculiar 
properties. Most of the numerical results presented in this paper concern model 
3. This model revealed to be much more stable under changes of its free 
parameters than the old GOY model. Why from different (severe) truncations 
of Navier-Stokes eqs. one ends up with so different dynamical behaviors is 
certainly the most stimulating question arising from our study. 
As for model 2, our preliminary numerical simulations suggest the presence of 
a relevant backward energy transfer leading to possibly strong deviations from 
the Kolmogorov scaling. A detailed study of model 2 will be reported later 

The crucial role played by inviscid invariants seems to be confirmed, specially 
when the phenomenological analysis suggests the presence of two simultaneous 
transferred quantities. 

Let us notice that usual arguments based on the analysis of absolute equilib- 
rium behavior [^, |2T| in order to extrapolate strongly dissipative effects like 
intermittency seems to fail in this class of shell models. Indeed, starting from 
the analysis of absolute equilibrium one should conclude that by changing a, 
and therefore the dimension of the second invariant, also the inertial properties 
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in the dissipative case should change. This is definitely true in the GOY model 
but definitely false in model 3. 

We have also written down the PDEs describing the semi-universal (equal for 
models 1,2 and 3) continuum limit. Investigations of this PDE are in progress. 
Let us conclude with some speculation. Such a rich behavior of these four 
models naturally suggests that could be important to ask the same question in 
the original Navier-Stokes fiow: what are the characteristics of the four sub- 
models obtained in the full Navier-Stokes decomposition taking into account 
separately the four different helicity classes? Can one recognize also in the case 
of full 3D dynamics sub-models with important back transfer of energy (like 
shell model 4 and shell model 2), or sub-models with different sensibility to the 
presence of helicity (like shell model 1 and shell model 3)? 
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FIGURE CAPTIONS 



Figure (1): Energy exchange in the one-triad system for the four models. 
Dashed (solid) arrows point towards the mode that receives less (more) 
energy. 

Figure (2): Variations with a of the energy rates in the one-triad system 
for the four models. For each model are shown the E of the two modes 
that receive energy from the unstable one, whose energy rate is always 
kept equal to 1. 

a): E2 (solid line) and Ej, (dashed line) vs a for model 1 (E\ = 1); b): 
El (solid line) and E^ (dashed line) vs a for model 2 {E2 = 1); c): E2 
(solid line) and E3 (dashed line) vs a for model 3 {Ei = 1); d): Ei (solid 
line) and E^ (dashed line) vs a for model 4 {E2 = 1). 

Figure (3): Helicity exchange in the one-triad system for the four models. 
Dashed (solid) arrows point towards the mode that receives less (more) 
helicity. 

Figure (4): The logarithm log2(S'p(n)) of the structures functions of model 
3 vs \og2{kn). The parameters values are a = 1 and A = 2. 

Figure (5): The C{p)s of GOY model (from |[T^) and model 3. The 
parameters values are a = 1 and A = 2. Error bars for data concerning 
model 3 take into account both statiscal and power-law fit errors. 

Figure (6): The C{p)s of GOY model and model 3 for different parameters 
sets (a, A). A is always kept equal to 2. The C(p)s (p = 1, ••, 7) for the 
GOY case (1,2) are taken from ||l^. Notice that for model 3 all data 



sets collapse on the same curve for different a values. Error bars for data 
concerning model 3 take into account both statiscal and power-law fit 
errors. 

Figure (7): The ({p)s of GOY model and model 3 for different parameters 
sets (a. A), a is always kept equal to 1. The C(p)s {p = 1, ••, 7) for the 
three GOY cases are taken from [|T^. Notice that for one value of A 
(A = 1.5) our numerically evaluated C{p)s are slightly different from those 
found in the GOY model and in model 3 for other A values. We are not 
able to conclude whether this discrepancy is important or not. Error bars 
take into account both statiscal and power-law fit errors. 



20 







S2 


S3 


S4 


S5 


S6 


1 


+ 










+ 


2 






+ 




+ 




3 




+ 




+ 






4 


+ 


+ 


+ 


+ 


+ 


+ 



Table 1: Helicity indeces of equations (|T^) for the four models. 
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Table 2: Coefficients of equations (18) for the four models 
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Table 3: Factors in the equations for the four models. 
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